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Abstract 

We consider the evolution of logistic maps under long-term memory. The memory effects are 
characterized by one parameter a. If it equals to zero, any memory is absent. This leads to the 
ordinary discrete dynamical systems. For a = 1 the memory becomes full, and each subsequent 
state of the corresponding discrete system accumulates all past states with the same weight just as 
the ordinary integral of first order does in the continuous space. The case with < a < 1 has the 
long-term memory effects. The characteristic features are also observed for the fractional integral 
depending on time, and the parameter a is equivalent to the order index of fractional integral. We 
study the evolution of the bifurcation diagram among a = and a = 0.15 . The main result of 
this work is that the long-term memory effects make difficulties for developing the chaos motion in 
such logistic maps. The parameter a resembles a governing parameter for the bifurcation diagram. 
For a > 0.15 the memory effects win over chaos. 

PACS numbers: 05.40.-a, 05.45.-a, 05.60.-k, 82.40.Bj 
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The treatment of nonlinear dynamics in terms of discrete maps (difference 
equations produced by numerical methods) is a very important step in study- 
ing the qualitative behaviour of continuous systems described by differential 
equations. The logistic map represents one of the most important examples of 
an one-dimensional discrete nonlinear map with the bifurcation scenario well 
known. The complicated behaviour exhibited by the logistic map is typical for 
a whole class of dynamical systems. Fractional calculus occupies an appreciable 
place in order to describe various kinds of wave propagation in complex media, 
fractional kinetics of Hamiltonian systems, anomalous diffusion and relaxation, 
etc. The fractional operator is a natural generalization of the ordinary differen- 
tiation and integration. When the operator depends on time, it is characterized 
by long-term memory effects. The effects correspond to intrinsic dissipative pro- 
cesses in the physical systems. In the application to discrete maps this means 
that their present state evolution depends on all past states with a power weight. 
Appearance of long-term memory effects in the logistic map makes the coupling 
among states stronger. This feature is plainly directed against the development 
of the chaotic dynamics. 



I. INTRODUCTION 



The evolution of discrete dynamical systems is described by the variables measured in 
discrete time steps. The behavior of such systems is governed by return maps. They connect 
the (n + l)th value of variables with the preceding nth value of the variables through a func- 
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tional dependence. A discrete dynamical system can demonstrate a complex development 
even in simple systems with one variable l|. If the system falls into a series of bifurcations, 
it exhibits a transition from a variety of periodic cycles to chaos. 

A statistical treatment of the macroscopic equation of motion leads to memory effects 
in dynamical Onsager coefficients, if Markovian-like approximations are not admissible 
The application of memory effects to the discrete systems signifies that their dynamics at 
time t + 1 will depend only not on time t, but on former times t — l,t — 2,t — 3 and so 
on. Here time is equivalent to the number of step. Therefore, n + 1 can be substituted 
for t + 1 as well as n can be t, and t — 1, t — 2, . . . can be replaced by n — 1, n — 2, . . . 
respectively. The introduction of a retardation of the linear term in the logistic equation 
modifies the Feigenbaum scenario |3j. In this case the periodic orbits are shifted. Several 
forms of nonlinear maps with memory were also considered by a formal way in J]. 

Last time the interest to the study of long-term memory effects and fractional kinetics in 



physical systems increased too much 



-|8j. The fractional operator with respect to time is 



characterized by long-term memory effects, whereas one with respect to coordinates possesses 



non-local (long-range) effects 



10|. The direct relationship among the long-term memory 



effects, fractional calcu 
been established in 11 



us and the stable distributions from the theory of probability has 



12[. Then it was shown 



13[ that perturbed by a periodic force, the 



nonlinear oscillator with fractional derivative exhibits a chaotic motion called the fractional 
chaotic attractor. Although the fractional nonlinear oscillator behaves like the stochastic 
attractor in phase space, being periodically perturbed, the role of the polynomial dissipation 
leads to a degradation of the attractor structure. The aim of this work is to present the 
analysis of the influence of long-term memory effects on the behavior of discrete systems. 
The paper starts with a discrete conception of long-term memory effects (Section HT|) . 



They are characterized by the parameter similar to the order index of fractional integral. By 
means of the computer simulation of bifurcation diagrams for the quadratic and triangular 
maps (Section IIIip with the memory contribution, we demonstrates how the memory effects 
strangle the chaotic motion of discrete dynamical systems. 

II. MATHEMATICAL DESCRIPTION OF LONG-TERM MEMORY 

The mapping x n+ i = f(x n ) does not have any memory, as the value x n+ \ only depends 
on x n . The introduction of memory means that the discrete value x n+ i is connected with 
the previous values x n , x n _i, x n _ 2 , . . ., X\. Particularly, any discrete system will have a full 
memory, if each state of the system is a simple sum of all previous states: 

n 

x n+l — f{ x i)i 
i=l 

where f(x) is an arbitrary function suitable for the discrete map. However, the above 
expression may tend to infinity because of the sum accumulating all the values. In order to 
have a fixed point in this mapping, the expression should be normalized as 

1 n 

x n+x = - J2 f( x i) ■ (!) 
When p is a fixed point, the map (JT]) gives 

1 - 

The starting conjection is that the map under long-term memory is expressed in terms of 

I - 

x n+l — — / y bjf(Xi), 

II i=l 

where the weights &j and the parameter a characterize the non-ideal memory effects. 
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Let the weights take the form 

' (1 + a)n a - n a+1 + (n - l) a+1 , if i = 0; 

1 , if i — n; 

{n-i + l) a+1 -2{n-i) a+l + 

t +(n - i - , ifO<i<n-l. 
In this connection it should be pointed out that the similar representation is used for the 

robast algorithm of numerical fractional integration 14]. It has found an interesting appli- 



cation for the analysis of the nonlinear reaction with fractional dynamics 



the study of fractional oscillations 



15| as well as for 



161 ] . Then the normalized mapping is written as 



x n+l 



(2) 



Here the sum starts with i = (though it could be realized from i — 1) by reason of the 
simplicity in writing the weights cj . To calculate the normalization factor presents no 
difficulty, namely 



i=0 



(n) 



1 + O-V 



It turns out that the factor is a power function just as one in our foregoing conjection. 

Now consider the following limit cases. If a = 0, the memory is absent. Really, in this 

case all the weights are equal to zero with the exception of = 1. Thus, we come to an 

ordinary map x n+ i = f(x n ). For a = 1 the memory effects will be full, namely 

fl, if i = 0; 

1 , if i = n; 

v 2, if < i < n — 1. 
Then this map becomes x n+ i = ^[(/(a^o) + f(x n ))/2 + J27=i f( x i)]- The full memory is ideal 
because each new state of the discrete system sticks to the system's memory and has the 
same action upon next states as all the others in memory. When < a < 1, the memory 
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effects have a long-term dependence so that 

f n 

^ +1 = irT ^gc i ( "V(^). (3) 

Certainly the memory is not ideal, as the contribution of more early states is noticeably 
less than the contribution of following ones on the present state of such systems. This is 
something like a part of system states being lost. Recall that if a system does not remember 
any previous state except for the present, it has no memory effects. 

In the next section we will show the results of numerical simulations for the triangular 
and quadratic maps under the long-term memory effects. This allows us to estimate their 
influence on the evolution of the discrete systems. The memory effects are governed by 
means of the parameter a, whereas the parameter of type r serves as an order parameter 
for the onset of chaos. 



III. NUMERICAL RESULTS 



The triangular map is expressed in terms of the function 



A(x) = r(l - 2 



x 



and the quadratic map is written as 



where r, q are constants. These logistic maps are very popular in the theory of deterministic 
chaos. They help easily to understand main features of discrete dynamical systems. In 
particular, the quadratic map is the simplest nonlinear difference equation, appears in many 
contexts, for example, a strongly damped kicked rotator or the growth model of a population 



in a closed area 



1]. For the triangular map with r < 1/2 the origin x = is the only stable 
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FIG. 1: The bifurcation diagrams for the triangular logistic map with long-term memory effects. 
The parameter a accounts for the memory contribution, the constant r relates to the onset of chaos 
for this map (all details in the text). 

fixed point to which all points in the interval [0 1] are attracted. The value r plays a role 
of "order parameter", which indicates the onset of chaos. For r > 1/2 two unstable fixed 
points emerge. If the magnitude of the parameter r increases still more, then the information 
about the position of a point in [0 1] is lost. This may result in chaos. The function A(x) is 
a simple model which for r > 1/2 generates chaotic sequences x , A(x ), A[A(x )}, . . ., and 
due to its simple form, all quantities in this chaotic state can be calculated explicitly. We 
go the same way for the maps taking into account the long-term memory effects. 

The long-term memory distinguishes itself by a slow decay contribution of former values 
for discrete maps. The procedure of calculating the string of system states requires to 
evaluate in the coefficients . Therefore, the duration of calculations notably increases 
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FIG. 2: The bifurcation diagrams for the quadratic logistic map with long-term memory effects. 
The parameter a characterizes the memory influence, the value q is the constant of the map (see 
details in the text). 

with the growth of n. To remember "the simple dynamical systems do not necessarily lead 



to simple dynamical behavior" [17|], we guess that a weak memory effects with a — > 
almost will not influence on the complicated behavior of such discrete systems. By means of 
numerical simulation we establish the evolution of the simple logistic maps under long-term 
memory, when the parameter a tends to one. For this purpose the bifurcation diagram is 
constructed. Figured] shows a set of bifurcation diagrams to the triangular map with various 
values a. A number of diagrams, represented in Fig. [2j corresponds to the quadratic case. 
The number of iterations n equals to 3500. In these figures it is seen that the memory effects 
step-by-step in the growth of a block up the development of chaotic behavior in the logistic 
maps. After a > 0.155 the evolution of the given maps tends to a regular stable fixed point. 
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What can we say about a sequence of iterates if there are unstable fixed points? The nth 
iterate A n (x) is piecewise linear and has the slope -^A n (x) = 2 n , except for the countable 
set of points j ■ 2~ n with j = 0, 1, ... , 2 n . The Liapunov exponent measures the average 
loss of information about the position of a point in [0 1] after one iteration. For the general 
triangular map under long-term memory effects the Liapunov exponent becomes 

X(x) = lim -l n y Cl (n) -^A l (a;) < lim - In V cf f A ! (x) 

in in 

= «& " <=!"'(»■)'} = to - ln{(2r)» £ c<"»(2r)'-} 

n i=0 " i=0 

= ln(2r). 

The equality A(x) = ln(2r) holds true for a = 0, when any memory effects are missing. If 
< a < 1, then the Liapunov exponent decreases. Therefore, the chaotic states typical 
for the ordinary triangular map die out with the growth of memory effects, i. e. when the 
parameter a tends to one. 

IV. CONCLUSIONS 

We have shown that if a discrete dynamical system is exposed to the long-term memory 
effects, then the chaotic sequences which are presented in the system without any memory 
are squeezed out. As the memory mounts, chaos and its traces disappear. The magnitude 
a crit m 0.155 looks like a "critical point" that indicates a transition from chaos to order. 
It should be noticed that the fine structure in the iterates for the maps is just washed out 
because of the memory effects. Nevertheless, for sufficiently small values of a there is still a 
remarkable transition to chaos. 
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